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Editor’s note 


It is great that Treubia volume 37 can be published in year 2010. 
Recently, it was difficult to get appropriate papers since animal taxonomy has 
not been an attractive subject in the field of biology. There was a lack of 
submitted manuscripts in 2009 that made Treubia could not be published in 
year 2009. 

This volume of TREUBIA contains five papers of vertebrates and 
invertebrates. Three papers (nematode, rats and land snail) were from the results 
of field works in eastern part of Indonesia i.e. West Papua which was rarely 
explored. 

Also, this year Indonesian zoologist’ community lost the pioneer and 
expert in parasite taxonomy, Dr. Sampurno Kadarsan. His name has been used 
to name new species of leeches, tick, rat, lizard and frog by his successors to 
acknowledge his impact and contribution. He served as an editor of Treubia 
from 1992 to 1997 and was a proof reader for some years until his permanent 
retarded eye sight. So, his death was a great lost for all of us especially for the 
Museum Zoologicum Bogoriense. 

Finally, I would like to thank all of the co-editors, referees, computing 
assistant, secretary and administrative assistant for their collaborative work. I 
acknowledge financial support from the Director of Research Centre for Biology 


LIPI to publish this precious journal. 


Cibinong, 15 December 2010 


Dewi M. Prawiradilaga 
Chief Editor 
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ABSTRACT 


The morphological study was focused on the analysis of morphometric 
variations of the spiny rats (Maxomys Jentink) across their Indonesian and Malaysian 
ranges in order to clarify the taxonomic status of this genus. The genus was widespread 
throughout South Asia, and many taxa (species and subspecies) have been described. 
Univariate and multivariate analyses (Principal Component Analyses and Discriminant 
Function Analyses) were carried out based on 24 craniometric characters in 187 skulls. 
Results indicated that Maxomys consist of 13 species: M. whiteheadi, M. surifer, M. 
rajah, M. ochraceiventer, M. baeodon, M. alticola, M. musschenbroekii, M. hellwaldii, 
M. wattsi, M. pagensis, M. bartelsii, M. inas and one distinct population of M. 
whiteheadi from 16 known species of Maxomys in Malaysia and Indonesian 
archipelago. The multiple regression analyses showed no sexual dimorphism in relation 
to the characters used; however all characters were geographically significant in some 
species. 


Key words: geographical variation, Maxomys, morphology, statistical analyses, 
taxonomy. 
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from Indonesia and Malaysia based on morphological study 


INTRODUCTION 


Recent reappraisal of the mammalian fauna of the Indo-Malayan region 
has pointed to a dramatically different view of the biogeographic boundaries 
of the Indo-Malayan and Australian regions from the traditional ones based on 
Wallace’s and Weber’s lines (Kitchener et al. 1990. Kitchener & Maryanto 
1995). The Indo-Malayan region presents a unique opportunity to study either 
intra or inter specific morphological variations in relation to geographical 
barriers. Comprehensive and accurate re-evaluation of morphological characters 
in certain cases resulted in the discovery of new species. Such re-evalution 
has been conducted in Indonesia for bats, rodents and primates (Kitchener & 
Maryanto 1993, Maryanto et. al. 1997, Merker & Groves 2006). 

Spiny rats, well known as Maxomys, are common rats which can be 
found in lowland and mountain forests up to 3000 m a.s.l., from primary forests, 
secondary forests, plantations and disturbed areas adjacent to forests (Payne et 
al. 1985, Corbet & Hill 1992, Francis 2008). The 16 known species of Maxomys 
have widespread distributions starting from the mainland of Southeast Asia, 
throughout much of the Indochinese region, Malaya Peninsular to Borneo, 
Sulawesi, and Palawan, as well as on several smaller islands of the Sunda 
Shelf (Ruedas & Kirsch 1997, Corbet & Hill 1992, Musser & Carleton 2005). 

A number of previous studies have been conducted in order to reveal 
the taxonomic status of Maxomys. However, many scientists have suggested 
that detailed and accurate re-evalution are still needed before the systematics 
of this genus is stable. Ruedas & Kirsch (1997) noted that Maxomys is a species- 
rich genus embracing a wide range of morphological, geographic, ecological, 


and altitudinal variations, and that most species still require rigorous definition. 
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Musser & Newcomb (1983), Corbet & Hill (1992), Watts & Baverstock 
(1994), Musser & Carleton (2005) all showed that Maxomys was not closely 
related to Rattus as earlier authors (e.g. Ellerman, 1941) had placed it. Musser 
et al. (1979), Corbet & Hill (1992), Musser & Carleton (2005) made detailed 
definitions of the Sundaic species of Maxomys, and offered a working hypothesis 
about the taxonomy and geographical distribution of Maxomys. The genus is 
divided into 16 species: M. alticola, M. bartelsii, M. baeodon, M. dollmani, M. 
hellwaldii, M. hylomyoides, M. inas, M. inflatus, M. moi, M. musschenbroekii, 
M. ochraceiventer, M. pagensis, M. panglima, M. rajah, M. surifer and M. 
wattsi. There was no systematic revision of this genus since then. This re- 
evaluation of the morphological characters of Maxomys was conducted with 
the main goal to reveal the morphological and geographical variations within 
this genus. The hypothesis tested here is whether there is significant variation 


in characters within Maxomys species from Indonesia and Malaysia 


MATERIALS AND METHODS 


A total of 187 skulls from adult specimens consisting of M. whiteheadi, 
M. rajah, M. ochraceiventer, M. surifer, M. baeodon, M. pagensis, M. bartelsii, 
M. inas, M. musschenbroekii, M. hellwaldii, M. wattsi, and M. alticola were 
studied (Appendix 1). All specimens studied were deposited at Museum 
Zoologicum Bogoriense (MZB) Indonesia, Museum Zoological UNIMAS 
(MZU) Malaysia and Raffles Museum, Singapore. Specimens were judged as 
adult based on fusion of basisphenoid or basioccipital sutures. Twenty four 
measurements of craniometric characters were made using digital calipers 


(Mitutoyo™) to 0.01 mm. The dimensions and limits are illustrated in Figure 1 
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and defined in Maryanto (2003). The characters were: GSL (greatest skull 
length), POW (post orbital width), ZB (zygomatic width), BoZP (breadth of 
zygomatic plate), IOB (interorbital breadth), NL (nasal length), NW (nasal 
width), BB (breadth of braincase), BH (braincase height), DL (length of 
diasthema), PL (palatal length), UMR (upper molar tooth row), IFL (length of 
incisive foramina), IFB (breadth of incisive foramina), M'W (upper molar 1 
width), M?W (upper molar 2 width), M*W (upper molar 3 width), М'М! 
(distance of right-left upper molar 1), M’M? (distance of right-left upper molar 
2), ММЗ? (distance of right-left upper molar 3), BL (bulla length), MSW 





Figure 1. Views of cranium and molars of Maxomys, showing limits of cranial and 
dental measurements. See text for abbreviations and additional informa- 
tion 
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(mesopterygoid fossa width), RAP (ramus angular process), MaL (length of 
mandible). 

Two way statistical analyses were used, univariate and multivariate 
analyses, as described in Maryanto & Sinaga (1998) and Maryanto (2003). 
Principal Component Analyses (PCA) and Discriminant Function Analyses 
(DFA) were initially run for all characters for all species as separated groups. 
Then DFA were run for a subset of five or more characters, based on criteria of 
minimised Wilk’s lambda to select the best discriminating variables. The results 
obtained with a reduced character set were very similar to those based on all 


characters (Kitchener et al. 1993). 


RESULTS AND DISCUSSION 


Univariate Statistics 

Mean, standard deviation, minimum and maximum values and samples 
size were represented in Table 1. These suggest that M. hellwaldii is the 
largest in size for all skull characters and M. whiteheadi (Small size or S) is 
the smallest species for most characters except for ZP, UMR, IFL, M1 W, M?W, 
M?W and BL larger than M. baeodon and NW, HB, D, and MSW slightly 
larger than M. inas. The means of M. whiteheadi Large size (L) are slightly 
smaller than M. ochraceiventer and clearly larger than M. whiteheadi Small 


size (S) and M. baeodon. 
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Multiple Regression 

Multiple regression analysis was employed to determine the influence 
of sex in variation. All of the characters examined were not significantly 
(P<0.05) influenced by sex (Table 2). Consequently, multivariate analyses were 


combined for males and females. 


Table 2. Multiple regression on sex, locality and interactions among these factors in 
Maxomys for skull and dentary characters. F values are presented for main 
effects with significance levels are *) 0.05>p>0.01; **) 0.01>p>0.001; and 
ээж) P20.001 


Dependent 





Variable Sex Locality Interaction 
GSL 1.297 12.062*** 0.260 
BB 0.099 9.372*** 0.414 
ZB 0.504 12.616*** 0.365 
ZP 1.711 8.983 *** 0.971 
IOB 1.540 16.442*** 0.430 
NL 0.433 13.686*** 0.302 
NW 0.964 11.050*** 0.107 
POW 0.770 0.858 *** 0.592 
HB 0.679 11.679*** 0.586 
р 0.523 8.980*** 0.158 
PL 0.955 7.474*** 0.183 
UMR 1.878 11.147*** 1.141 
IFL 0.547 15.472*** 0.320 
IFB 2.291 4.879*** 0.634 
MIW 1.398 12.826*** 0.449 
M2W 0.827 9.988 *** 0.541 
M3W 0.224 4.17 | *** 0.873 
MIMI 0.157 6.658 *** 0.502 
M2M2 0.336 6.775*** 0.814 
M3M3 1.150 11.148*** 0.697 
BL 0.036 14.501 *** 0.613 
MSW 0.027 13.255*** 0.899 
RAP 0.725 10.766*** 0.266 
DL 0.362 12.986*** 0.497 
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Multivariate Analyses 


Morphological Variation Within Genus Maxomys 

PCA analysis of morphology Maxomys indicates significant variation 
between or within M. whiteheadi (L), M. whiteheadi (S), M. rajah, M. surifer, 
M. baeodon, M. ochraceiventer, M. alticola, M. musschenbroekii, M. hellwaldii, 
M. wattsi, M. pagensis, M. bartelsii and M. inas. Five principal components 
(PCs) extracted from 24 measurements accounted for 90,0 % of the total 
variation (Table 3). The DFA is based on 24 characters, and repeated using a 
reduced set of five characters in order to avoid over fitting the data (Kitchener 
& Maryanto 1995). The five characters were selected by minimising Wilk's 
lambda on the first canonical variate. The DFA plot is based on five characters 
similar to the complete characters set (Maryanto 2003). The five characters in 


the DFA are GSL, BB, ZB, NL, PL (Table 4). 


Table 3. Total variation accounted for 24 measurements on five principal components 


(PCs) 
Initial Eigenvalues шилэн Squared 
% of Cumulative % of Cumulative 
Total variance % Total variance % 
1 1771 73.80 73.80 17.71 73.80 73.80 
2 144 6.01 79.81 1.44 6.01 79.81 
3 1.08 4.52 84.32 1.08 4.52 84.32 
4 082 3.41 87.73 0.82 3.41 87.73 
5 054 2.27 90.00 0.54 2.27 90.00 


Extraction Method: Principal Component Analysis. 
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Table 4. Standardised and unstandardised (in the brackets) canonical variate function 
coefficients derived from analysis of five characters 


1 2 
GSL 0.017(0.013) — 0.830(0.653) 
BB 0.461(0.984) -0.108(-0.230) 
ZB 0.005(0.010) 0.091(0.164) 
М, 0.477(0.754) 0.433(0.686) 
PL 0.476(0.886) -1.212(-2.254) 
(Constant) (-36.865) (0.182) 


Function 


3 4 
-0.603(-0.474) 1.100(0.866) 
-0.496(-1.059) 0.423(0.904) 

1.436(2.581) -0.176(-0.316) 
-0.196(-0.311) -0.963(-1.524) 
-0.048(-0.090) -0.290(-0.539) 


(-3.679) (-12.610) 


5 


-0.794(-0.625) 
0.762(1.626) 
0.253(0.454) 
0.149(0.236) 

-0.236(-0.439) 


(-5.679) 


This analysis extracted five significant functions, which accounted 


for 10096 of total variations. Discriminant Functions (DF) 1, 2, 3, 4, and 5 
accounted for 88.49 96 (P«0.001; df : 60), 4.92 96 (P«0.001; df : 44), 3.23 96 
(P«0.001; df : 30), 1.88 96 (P«0.001; ағ: 18) and 1.49 96 (P«0.001; df : 8) of 


the variation between populations respectively (Table 5) with 87.7 96 of the 


originally grouped cases correctly classified and 12.346 misclassified. 


Furthermore, the number of functions at group centroids between skull, dental 


and dentary characters from each species are distinct. 


Table 5. Total explained variation for 24 measurements on five Discriminant functions 


A шо N н 


(DFs) 
Function Eigenvalue % of Cumulative Canonical 
Variance % Correlation 
21.38(a) 88.49 88.49 0.98 
1.19(a) 4.92 93.41 0.74 
0.78(a) 3.23 96.63 0.66 
0.45(a) 1.88 98.51 0.56 
0.36(a) 1.49 100.00 0.51 


a. First 5 canonical discriminant functions were used in the analysis. 
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Figure 2. 


Figure 3. 
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Factor 1 








Factor 2 
Plot of scores on factor axes 2 and 1 from principal component analysis 
(PCA) of 24 skull measurements of Maxomys. Groups classification: a, 
M. whiteheadi (S); b, M. surifer; c, M. whiteheadi (L); d, M. rajah; e, M. 
ochraceiventer; f, M. baeodon; g, M. alticola; h, M. musschenbroekii; i, 
M. hellwaldii; |, M. wattsi; k, M. pagensis; 1, M. bartelsii; m, M. таз; (1), 
large body size; (2), medium body size; and (3), small body size 





Factor 1 





Factor 3 


Plot of scores on factor axes 3 and 1 from principal component analysis 
(PCA) of 24 skull measurements of Maxomys. Groups classification: a, 
M. whiteheadi (S); b, M. surifer; c, M. whiteheadi (L); d, M. rajah; e, M. 
ochraceiventer; f, M. baeodon; g, M. alticola; В, M. musschenbroekii; i, 
M. hellwaldii; |, M. wattsi; k, M. pagensis; 1, M. bartelsii; m, M. таз; (1), 
large body size; (2), medium body size; and (3), small body size 
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Factor 1 








Factor 4 


Figure 4. Plot of scores on factor axes 4 and 1 from principal component analysis 
(PCA) of 24 skull measurements of Maxomys. Groups classification: a, 
M. whiteheadi (S); b, M. surifer; c, M. whiteheadi (L); d, M. rajah; e, M. 
ochraceiventer; f, M. baeodon; g, M. alticola; h, M. musschenbroekii; i, 
M. hellwaldii; j, M. wattsi; k, M. pagensis; 1, M. bartelsii; m, M. inas; (1), 
large Боду size: (2). medium bodv size: and (3). small body size 





Factor 1 








Factor 5 


Figure 5. Plot of scores on factor axes 5 and 1 from principal component analysis 
(PCA) of 24 skull measurements of Maxomys. Groups classification: a, 
M. whiteheadi (S); b, M. surifer; c, M. whiteheadi (L); d, M. rajah; e, M. 
ochraceiventer; f, M. baeodon; g, M. alticola; h, M. musschenbroekii; i, 
M. hellwaldii; |, M. wattsi; k, M. pagensis; 1, M. bartelsii; m, M. таз; (1), 
large body size; (2), medium body size; and (3), small body size 


66 


Treubia 2010, 37: 49 - 82 


Function 1 





Function 2 


Figure 6. Plot of scores on function axes 2 and 1 from discriminant function of five 
dominant variables skull measurements of Maxomys. Groups classification: 
a, M. whiteheadi (S); b, M. surifer; c, M. whiteheadi (L); d, M. rajah; e, M. 
ochraceiventer; f, M. baeodon; g, M. alticola; h, M. musschenbroekii; i, M. 
hellwaldii; j, M. wattsi; k, M. pagensis; 1, M. bartelsii; m, M. inas; (1), 
large body size; (2), medium body size; and (3), small body size 


Function 1 











Function 3 


Figure 7. Plot of scores on function axes 3 and 1 from discriminant function of five 
dominant variables skull measurements of Maxomys. Groups classification: 
a, M. whiteheadi (S); b, M. surifer; c, M. whiteheadi (L); d, M. rajah; e, 
M. ochraceiventer; f, M. baeodon; g, M. alticola; h, M. musschenbroekii; 
i, M. hellwaldii; |, M. wattsi; К, M. pagensis; 1, M. bartelsii; m, M. inas; 
(1), large body size; (2), medium body size; and (3), small body size 
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Function 1 











Function 4 


Figure 8. Plot of scores on function axes 4 and | from discriminant function of five 
dominant variables skull measurements of Maxomys. Groups classification: 
a, M. whiteheadi (S); b, M. surifer; c, M. whiteheadi (L); d, M. rajah; e, M. 
ochraceiventer; f, M. baeodon; g, M. alticola; h, M. musschenbroekii; i, M. 
hellwaldii; j, M. wattsi; К, M. pagensis; 1, M. bartelsii; m, M. inas; (1), 
large body size; (2), medium body size; and (3), small body size 
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Function 5 


Figure 9. Plot of scores on function axes 5 and 1 from discriminant function of five 
dominant variables skull measurements of Maxomys. Groups classification: 
a, M. whiteheadi (S); b, M. surifer; c, M. whiteheadi (L); d, M. rajah; e, M. 
ochraceiventer; f, M. baeodon; g, M. alticola; h, M. musschenbroekii; i, M. 
hellwaldii; |, M. wattsi; К, M. pagensis; 1, M. bartelsii; m, M. inas; (1), 
large body size; (2), medium body size; and (3), small body size 
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The configurations of PCA and DFA plots suggest that spiny rats in 
Indonesia and Malaysia comprise thirteen major morphological distinctness 
(figures 2; 3; 4; 5; 6; 7; 8; 9). These are: M. rajah, M. surifer, M. baeodon, M. 
alticola, M. ochraceiventer, M. whiteheadi (S), M. musschenbroekii, M. 
hellwaldii, M. wattsi, M. pagensis, M. bartelsii, M. inas and the two distinct 
populations of М. whiteheadi (L). Plots of PCA and DFA show that the genus 
Maxomys is divided into three major groups based on body size: group 1 has 
large body size and contains four species: M. rajah, M. pagensis, M. wattsi, 
and M. hellwaldii. Group 2 contains medium body sized species: M. surifer, 
М. alticola, M. bartelsii, М. musschenbroekii, M. ochraceiventer and M. 
whiteheadi (L). Group 3 has small body size and comprises M. whiteheadi 
(S), M. baeodon and M. inas. 

Maxomys whiteheadi (L) has distinct morphological characters from 
the other Maxomys shown by good separation from other species and large 
variation in 24 of the skull, dental and dentary characters measured. This species 
seems similar in morphological features to M. whiteheadi (S) and M. 
ochraceiventer, supported by overlapping measurements in certain individuals 
of skull characters. Characters examined of this species are strictly distinct 


from M. ochraceiventer collected from entire Borneo island. 


Geographical variation within Maxomys whiteheadi Small and Large based 
on morphological characters 

DFA analysis of geographical variation within M. whiteheadi (S) 
indicates significant variation from three different regions: Sumatra, southern 
region of Borneo, and northern region of Borneo. Two discriminant functions 


extracted accounted for a 100% of the total variation (Table 6). The DFA was 
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done on all 24 characters and repeated using a reduced set of seven characters 
in order to avoid over fitting the data (Kitchener & Maryanto 1995). The seven 
characters were selected by minimising Wilk’s lambda on the first canonical 
variate. The DFA plot is based on seven characters similar to the complete 
characters set (Maryanto 2003). The seven characters in the DFA are BB, 
ТОВ, NL, NW, M'M!, М?М?, and RAP (Table 7). This analysis extracted two 
significant functions, which accounted for 100% of total variations. 
Discriminant Functions (DF) 1 and 2 accounted for 94.0496 (P«0.001; df : 14) 
and 5.96% (P«0.05; df : 6) of the variation between populations respectively 
(Table 6) with 88.246 originally grouped cases correctly classified and 12.996 
misclassified; 82.9% of cross-validated grouped cases were correctly classified. 
Furthermore, a number of functions at group centroids between skull, dental 


and dentary characters from each species are also distinct. 


Table 6. Total explained variation for 24 measurements on two Discriminant functions 


(DFs) 
Function Eigenvalue % of variance Cumulative 96 enon! 
Correlation 
1 3.35 (a) 94.04 94.04 0.88 
2 0.21(a) 5.96 100.00 0.42 


Table 7. Standardised and unstandardised (in the brackets) canonical variate function 
coefficients derived from analysis of seven characters 


Function 
1 2 

BB 1.034(2.558) 0.059(0.145) 

IOB -0.691(-2.522) -0.462(-1.684) 
NL 0.820(1.413) 0.119(0.205) 

NW 0.334(1.663) 0.143(0.714) 

М1М1 -0.390(-2.044) 0.551(2.890) 

M3M3 -0.509(-2.000) 0.580(2.282) 

КАР 0.009(0.019) -0.944(-2.089) 
(Constant) (-28.725) (1.922) 


70 


Treubia 2010, 37: 49 - 82 


Function 1 





Function 2 


Figure 10. Plot of scores on function axes 2 and 1 from discriminant function of 
seven dominant variables skull measurements from M. whiteheadi (S). 
Group classifications: a, Southern Borneo; b, Sumatra; and c, Northern 
Borneo 


The configuration of DFA plots suggests that M. whiteheadi (S) is 
divided into three major geographical populations (Figure 10). The first group 
comprises populations from (1) Kutai, East Kalimantan and DAS Sebangau, 
Central Kalimantan; (2) Kayan Mentarang; North-eastern part of Kalimantan; 
and (3) Mount Palung, Ketapang, West Kalimantan. The second group 
comprises populations from Sarawak (Northern Borneo; Kubah National Park 
and UNIMAS peatswamp forest), and the third group comprises populations 
from the northern part (Deli and South Tapanuli) and the southeastern part 
(Kalianda, Lampung) of Sumatra. The populations from Sumatra are slightly 
similar in skull features and some characters overlap for some individuals from 


Northern or Southern Borneo. However, the populations from Northern Borneo 
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are Clearly distinct from Southern Borneo populations. This distinction possibly 
reflects geographical barriers such as mountains, highlands or large rivers in 
Borneo. 

The DFA plots for M. whiteheadi (L) show that there is geographical 
variation among the populations. Within this populations, there is significant 
variation between Borneo mainland (represented by specimens from Melak, 
Kalimantan, Sarawak, and Sabah) and those from both Sumatra and Peninsular 
Malaysia (Figure 11-13). 

Four discriminant functions were extracted which accounted for 100% 
of the total variation (Table 8). Discriminant Functions (DF) 1, 2, 3 and 4 
accounted for 92.05% (P<0.001; df : 32), 6.04% (Р-0.001, АГ: 21), 1.44% 
(Р<0.01; df : 12), and 0.47% (P>0.05; df : 5) of the variation between 


populations respectively with 100% of the originally grouped cases correctly 














а 
& 
104 
ва ® 
За ® a ^ 
- 04 i ау 
c 
o 
o 
= 
Ц. 4104 
204 
— "ds 
erm 29 
CH 
-30 
+ - e -— € --— 
$ 4 -2 0 2 4 6 


Function 2 
Figure 11. Plot of scores on function axes 2 and 1 from discriminant function of eight 
dominant variables skull measurements from M. whiteheadi (L). Group 
classifications: a, Southern Borneo; b, Sumatra Island; c, Malaya Penin- 
sular; d, Northwestern Borneo; and e, Northeastern Borneo 
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Function 1 





Function 3 


Figure 12. Plot of scores on function axes 3 and 1 from discriminant function of eight 
dominant variables skull measurements from M. whiteheadi (L). Group 
classifications: a, Southern Borneo; b, Sumatra Island; c, Malaya 
Peninsular; d, Northwestern Borneo; and e, Northeastern Borneo 


Function 1 





Function 4 


Figure 13. Plot of scores on function axes 4 and 1 from discriminant function of eight 
dominant variables skull measurements from M. whiteheadi (L). Group 
classifications: a, Southern Borneo; b, Sumatra Island; c, Malaya 
Peninsular; d, Northwestern Borneo; and e, Northeastern Borneo 
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Table 8. Total explained variation for 24 measurements on four Discriminant functions 


(DFs) 
Function Eigenvalue % of Cumulative % noe! 
Variance Correlation 
1 191.54(a) 92.05 92.05 1.00 
2 12.57(а) 6.04 98.09 0.96 
3 3.01(a) 1.44 99.53 0.87 
4 0.97(a) 0.47 100.00 0.70 


a first 4 canonical discriminant functions were used in the analysis. 


Table 9. Standardised and unstandardised (in the brackets) canonical variate function 
coefficients derived from analysis of eight characters 


Function 
1 2 3 4 
BB 0.109(0.243) -1.681(-3.747) 0.816(1.819) 0.638(1.421) 
ZP 3.814(22.640) 0.006(0.034) — 0.818(4.858) -0.351(-2.083) 
IOB 0.068(0.271) 1.913(7.664) 0.273(1.095) 0.296(1.186) 
D -4.702(-8.325) -1.367(-2.420)  0.505(0.895) -0.342(-0.606) 
UMR 2.969(13.792)  -0.467(-2.171) 0.054(0.251) 0.477(2.216) 
ТЕГ, 3.728(10.776) 1.514(4.377) -1.129(-3.263) 0.398(1.151) 
MSW 2.477(8.892) 1.514(-0.404)  0.308(1.106) -1.045(-3.750) 
КАР -2.401(-5.631) -0.112(2.112) -1.336(-3.134) 0.086(0.201) 
(Constant) (-111.376) (2.196) (-17.022) (-26.249) 


classified and 90.5% cross-validated grouped cases correctly classified. The 
eight characters were selected by minimizing Wilk’s lambda on the first 
canonical variate in order to avoid over fitting data. The eight characters in the 
DFA are BB, ZP, IOB, D, UMR, IFL, MSW, and RAP (Table 9). 

The configuration of DFA plots suggest that M. whiteheadi (L) is 
divided into five major geographical populations (Figures 11, 12, 13); these 


are (1) populations from Melak, Kutai, East Kalimantan, Central Kalimantan 
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(Southern Borneo); (2) populations from Sumatra; (3) a Sarawak group 
(Northwestern Borneo); (4) a Malay peninsular group; and (5) a Sabah group 
(Northeastern Borneo). Each of the groups showed good separation between 
them. Populations from Sumatra seem similar and closely related to the Southern 
Borneo populations. The Sarawak and Malay peninsula groups are clearly 
different from the other populations and stand as a single group. These different 
populations are most likely a reflection of geographical barriers such as 
mountains, highlands, large rivers and seas in this region which have isolated 
the populations. 

Based on this study, we summarized that Maxomys can now be divided 
into 16 species with M. whiteheadi (L) as a distinct new species yet to be 
described. Based on the features of skins, skulls and body size configuration 
M. whiteheadi (L) seems more similar to M. whiteheadi (S) than to M. 
ochraceiventer or M. baeodon; both of which is sympatric. Furthermore, skull, 
dental and dentary characters of M. musschenbroekii and M. wattsi from 
Sulawesi show similar features to Maxomys populations from Kalimantan, in 
particular populations of M. alticola, M. ochraceiventer, M. surifer and M. 
whiteheadi (L) an observation supported by overlapping measurements in 
some cases. Maxomys wattsi appears to be closely related to M. rajah which 
is supported by overlaps in measurements. Musser et al. (1979) observed that 
M. musschenbroekii was a small rat, more like M. alticola, M. inas, and M. 
ochraceiventer in body size. It was also more similar to those species, and to 
M. hylomyoides and M. whiteheadi in features of skin, skull, and teeth than to 
either M. hellwaldii or M. dollmani. Ellerman (1941) indicated that M. 
musschenbroekii and M. whiteheadi (S) were conspecific but there is no 


evidence, at least from the study of skins and skulls, to support that view. 
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The biogeographic boundaries of the species of Maxomys in the 
Southeast Asia region, do not follow the Wallace’s Line and show no separation 
into eastern and western species groups. This is supported by strong evidence 
about similarity in morphological features from skull characters between eastern 
and western Maxomys species and consistently confirms the observations 


suggested by Ellerman (1941) and Musser et al. (1979). 
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Appendix 1. List of samples used in analysis 


Maxomys whiteheadi (Small body size) 

MZB 29090 Female, MZB 29092 Female, MZB 29073 Male, MZB 29093 
Male, MZB 29094 Male, MZB 29098 Male, MZB 29099 Male, MZB 29100 
Male, MZB 29102 Male, MZB 29103 Female, MZB 29104 Female, MZB 
23645 Female, MZB 23642 Female, MZB 23641 Male, MZB 29059 Female, 
MZB 29061 NR, MZB 29064 Female, MZB 29071 Female, MZB 29081 
Female, MZB 23607 Female, MZB 23601 Male, MZB 23651 Male, MZB 
23647 Male, MZB 23638 Male, MZB 15384 Male, MZB 26247 Male, MZB 
22220 Female, MZB 22212 NR, MZB 22215 Female, MZB 26251 Male, 
MZB 22217 Female, MZB 22222 Male, MZB 26253 Male, MZB 26234 Male, 
MZB 26246 Male, MZB 22210 Male, MZB 11647 Female, MZB 13338 
Female, MZB 18347 Male, MZB 18370 Male, MZB 18371 Male, MZB 24986 
Female, MZB 5661 Male, MZB 15383 Male, MZB 15569 Female, MZB 15571 
Female, MZB 20560 Male, MZB 20569 Female, MZB 20580 Female, MZB 
20567 Male, MZB 20582 Male, MZB 20559 Male, MZB 20,574 Female, 
MZB 20,571 Male, MZB 20,562 Female, MZB 20,568 Male, MZB 24944 
Female, MZB 24957 Female, MZB 2071 Female, MZB 24931 Female, MZB 
15385 Female, MZU TK152846 Female, MZU 819 Male, MZU TK156129 
NR, MZU 1844 Male, MZU TK152353 NR, MZU 326 Male, MZU 1052 
Female, MZU TK152362 NR, MZB МА014 NR, MZB 5679 Male, MZB 5673 
Male, MZB 5663 Male, MZB 24893 Male, MZB 24963 Male, MZB 24898 
Male, MZB 13031 Female, MZB 15382 Male, MZB 23648 Male, MZB 29065 
Female, MZB 29067 Female, MZB 29069 Female, MZB 22233 Female, MZB 
29074 Male, MZB 29078 Male, MZB 29082 Female, MZB 29083 Female, 
MZB 29084 Male, MZU TK152862 NR, MZB 29076 Female. 
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Maxomys whiteheadi (Large body size) 

MZB 29058 Female, MZB 29062 Female, MZB 29066 Female, MZB 29068 
Female, MZB 29070 Female, MZB 29077 Male, MZB 29080 Male, MZB 
18482 NR, MZB 29086 Female, MZB 28961 Male, MZB 28962 Female, MZB 
28960 Female, MZB 28963 Male, MZB 28967 Male, MZB 28968 Male, MZU 
TK153703 Male, MZU TK153717 Male, MZU TK152349 NR, MZU 
TK152861 NR, MZU 471 Male, MZU 1618 Female. 


Maxomys wattsi 
MZB 12159 Male, MZB 12160 Male, MZB 12161 Male, MZB 12163 
NR 


Maxomys surifer 

MZB 4961 Male, MZB 4962 Female, MZB 4969 Female, MZB 4971 Female, 
MZB 4974 Male, MZB 4976 NR, MZB 4977 Female, MZB 4978 Male, MZB 
4979 Male, MZB 4983 Female, MZB 4984 Male, MZB 4985 Male, MZB 
4987 Female. 


Maxomys rajah 

MZU TK152850 Male, MZU TK152852 Male, MZB 14756 Female, MZB 
14774 NR, MZB 14779 Male, MZB 14752 Male, MZB 22496 Male, MZB 
22621 Female, MZB 23661 Male, MZB 23662 Male, MZB 23663 Female, 
MZB 27436 Not recorded. 


Maxomys pagensis 


MZB 22355 Female, MZB 2885 Female, MZB 17565 Male. 


Maxomys ochraceiventer 


MZB 23606 Male, MZB 23971 Male, MZB 23970 Female. 
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Maxomys musschenbroekii 

MZB 26935 Female, MZB 26942 Female, MZB 26943 Male, MZB 26945 
Male, MZB 26946 Male, MZB 26950 Female, MZB 26951 Male, MZB 26953 
Male, MZB 26955 Male, MZB 26956 Male, MZB 26960 Female. 


Maxomys hellwaldii 
MZB 26964 Male, MZB 26971 Female, MZB 26982 Male, MZB 26988 
Male, MZB 26991 Male, MZB 27163 Female. 


Maxomys bartelsii 

MZB 2188 Male, MZB 9537 Male, MZB10135 Male, MZB 10872 Male, 
MZB 10875 Female, MZB 10876 Female, MZB 10877 Male, MZB 12430 
Female, MZB 12462 NR, MZB 12919 Male, MZB 22525 Male. 


Maxomys baeodon 

MZU TK152835 Male, Raffles Museum 3312 Male, Raffles Museum 3533 
Male, Raffles Museum 3543 Female, Raffles Museum 3571 Female, Raffles 
Museum 3224 Female. 


Maxomys alticola 

Raffles Museum 4213 Male, Raffles Museum 4203 Male, Raffles Museum 
4214 Male, Raffles Museum 4013 Male, Raffles Museum 4223 Female, Raffles 
Museum 4215 Male. 


Maxomys inas 


М7В 13432 Female. 
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